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ABSTRACT 

We report the detection of SiS rotational lines in high-vibrational states as well as SiO and SiC 2 lines in their 
ground vibrational state toward IRC+10216 during the Atacama Large Millimeter Array Cycle 0. The spatial 
distribution of these molecules shows compact emission for SiS and a more extended emission for SiO and 
SiC 2 , and also proves the existence of an increase in the SiC 2 emission at the outer shells of the circumstellar 
envelope. We analyze the excitation conditions of the vibrationally excited SiS using the population diagram 
technique, and we use a large velocity gradient model to compare with the observations. We found moderate 
discrepancies between the observations and the models that could be explained if SiS lines detected are optically 
thick. Additionally, the line profiles of the detected rotational lines in the high energy vibrational states show 
a decreasing linewidth with increasing energy levels. This may be evidence that these lines could be excited 
only in the inner shells, i.e., the densest and hottest, of the circumstellar envelope of IRC-1-10216. 

Subject headings: astrochemistry — circumstellar matter — line: identification — molecular processes — 
stars: AGB andpost-AGB — stars: individual (IRC-r 10216) 


1. INTRODUCTION 

Silicon is mostly locked in SiS, SiO and SiC 2 in the circum¬ 
stellar envelope (CSE) of the carbon-rich star IRC-rl0216, 
as evidenced observati o nally and predicted by mod¬ 
els (iQlofsson et al.l I1982I iLucas et ^ 119951: lAgundez et al.l 
l2012lk These molecules are efficiently formed in the gas 
phase, close to the stellar photosphere as a consequence of 
chem ical process es enabled under thermodynamical equilib¬ 
rium (lTsujilll973h . In the dust formation region (~5-20R^), 
the Si-bearing molecules are likely to condense onto the dust 
grains due to their highly refractory nature. The silicon con¬ 
tained in the dust grains can form molecules through grain- 
surface reactions. Also, the interaction of shocks produced by 
the pulsation of the star with the dust grains can extract cer¬ 
tain amounts of silicon from the grains and incorporate that 
sili con into the gas-phase to re act and form other species (see 
e.g. lCastro-Carrizo et al.ll2001h . Beyond this region, the abun¬ 
dances of Si-bearing molecules are expected to decrease up 
to the outermost shells of the envelope, where the interstellar 
ultraviolet (UV) radiation field dissociates all the remaining 
molecules. 

Previous interferometer observations showed the spatial 
distribution of these molecules in IRC-1-10216. The SiS J=5 
- 4, J=6 - 5, /=8 - 7, J=9 - 8 and J=\2 - 11 bright- 
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ness distributions display a quasi-circular shape with a di- 
ameter of ~20'' elongated along the nebul ar axis (RA.~20°, 
iBieging & Tafallal [19931: iLucas et al.lll995h . Recent observa¬ 
tions with the Combined Array for Research in Millimeter- 
wave Astronomy (CARMA) of the SiS /=14 - 13 y =0 and 
v=l lines have been reported bv lEonfria et al.l (120141) . where 
the v=0 line shows a circular and compact brightness distribu¬ 
tion of ~2" and displays maser emission. The v=l brightness 
distribution shows a compact source centered at the star posi¬ 
tion. 

SiO J=5 - 4 v=0 brightness distribution maps carried 
out with the Submi llimeter Array (SMA) were reported in 
ISchoier et al.l (l2006h . They show circular symmetry with a di¬ 
ameter of ~6" at the systemic velocity of the source, which is 
-26.5km s“^ (^e.g. ICernicharo et al.ll200Ql). The J=6 - 5 v=0 
brightness distribution reported in lEonfria et al.l (1201 4l) dis¬ 
plays a quasi-circular symmetry centered at the position of 
the star, with a diameter of ~3" elongated along the nebular 
direction (NE-SW). 

SiC 2 observations carried out with the Plateau de Bure In¬ 
terferometer (PdBI) and CARMA, show a brightness distri¬ 
bution composed of: (/) an elongated co mpact component lo- 
cated at the innermost regions of the CSE (lEonfria et al .120141) 
and, (//) a hollow s hell structure located at ~15" from the star 
(ILucas et al.lll995h . The for mation mechanism for th is outer 
component was suggested in ICernicharo et al.l (120101) . where 
the reaction between Si and C 2 H 2 yielding SiC 2 , could be re¬ 
sponsible for the SiC 2 enhancement in the outer envelope. 

In this work we present the Cycle 0 observations carried out 
with the Atacama Large Millimeter Array (ALMA) toward 
IRC -I- 10216. We detected emission of SiS /=15 - 14 lines of 
vibrationally excited states, from v=0 up to v=7, and tenta¬ 
tively of v=8, 9 and 10. 7=15-14 lines of different isotopo- 
logues are also detected: ^^SiS (v=0-5), ^^SiS (v=0-4), Si^^S 
(v=0-3), Si^^s (v=0-4), 29si33s (v=0) and (v=0). We 

also detected emission of SiO J=6 - 5 (v=0-2), ^^SiO J=6 
- 5 (v=0) and several lines of SiC 2 in the ground vibrational 
state. 
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Fig. 1.— From top to bottom: SiC 2 J- 106,5 v=0, ^^SiO /=6 - 5 v=0 and ^^SiS /=15 - 14 v=0. Three velocity offsets with respect to the Vlsr 
of the source are shown: +0, +6 and +12km s“^ The position of the star is shown as a white star. The first contour corresponds to 5 cr (see Table[T) and the 
rest correspond to 15, 25, 50, 75 and 90% of the maximum intensity, except for ^^SiS /= 15-14 where the first contour is at 25 cr level (a large-scale artificial 
modulation appears below that level). The synthetic beam is represented in the top right corner of each of the three maps series. The intensity scale in Jybeam“^ 
is shown at the right edge of the figure. The orientation is explicitly shown in the bottom right box. 


2. OBSERVATIONS 

The observations were carried out with ALMA0 between 
2012 April 8 and 23 during Cycle 0. IRC+10216 was ob¬ 
served in the frequency range 255.3 to 274.8 GHz (band 6) 
covered by four different setups with a bandwidth of ~5 GHz, 
a channel spacing of 0.49 MHz and an effective resolution of 
0.98 MHz. Detailed information of each setup is summarized 
in Table[T] The observations were performed using sixteen 
antennas covering baselines up to 402 m that allowed us to 
obtain an angular resolution of -0.6". The shortest baselines 
used were -20 m which allow us to recover structures with 
a size up to <12". Two runs of 72 minutes each were per¬ 
formed, of which 26 minutes correspond to correlations on 
source. Further details about calibration and imaging restora- 

^ This paper makes use of the following ALMA data: 
ADS/JAO.ALMA#2011.0.00229.S . ALMA is a partnership of ESO 
(representing its member states), NSF (USA) and NINS (Japan), to¬ 
gether with NRC (Canada), NSC and ASIAA (Taiwan), and KASI 
(Republic of Korea), in cooperation with the Republic of Chile. 
The Joint ALMA Observatory is operated by ESO, AUI/NRAO and 
NAOJ. 


tion can be found in ICernicharo et al.l (l2013h . The coverage 
of the uv plane achieved with the setup 6 provides low con¬ 
tributions of the sidelobes (<10% of the primary beam) to the 
dirty beam. For the rest of the setups the uv coverage is worse 
and large contributions of the sidelobes (up to 20-30% of the 
primary beam) appear in the dirty beam. 

The continuum comes from a point-like source located 
at q;= 9‘’47”57?446 and (5=13°16'43;'86(J2000), which is in 
good agreement with the position of IRC-hl0216 measured 
with the Very Larg e Array (VLA) with 40 mas resolution 
(iMenten et al.ll2012h . We measured an intensity peak of 650 
mJy beam“^ with an uncertainty of -8 %. 

The calibration of the data was performed using CASAE3 
and data analysis with GILDASE3- 

CASA (Common Astronomy Software Applications) is 
a comprehensive software package used to calibrate, image, 
and analyze radioastronomical data from interferometers. See: 
http://casa.nrao.edu 

GILDAS is a worldwide software mainly used to process, re¬ 
duce, and analyze astronomical single-dish and interferometric ob¬ 
servations. It is maintained by the Institut de Radioastronomie Mil- 
limetrique (IRAM). See: http://www. iram. fr/IRAMFR/GILDAS 
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Fig. 2.— ^^SiO 7=6 - 5 v=0 velocity interval averaged maps. The velocity averaged interval is shown at the top left corner of each box. These velocities are 
offset velocities with respect to the Vlsr of the source (-26.5 km s“^). The position of the star is shown as a white star. The first contour corresponds to 5 cr (see 
Table[TJ and the rest correspond to 15, 25, 50, 75 and 90% of the maximum intensity. The intensity scale in Jybeam“^ is shown at the right edge of the figure. 
The synthetic beam is represented in the top right corner of the bottom right box. The orientation is explicitly shown in the bottom right box. 


TABLE 1 

Observational parameters. 


Setup # Frequency range RMS 
(GHz) (mJy beam 

3 269.9-274.8 6 

4 265.0-269.9 10 

5 260.2-265.0 17 

6 255.3-260.2 6 


Synthetic beam FOV’*' 

) (^^x^O _ (^0 

0.61x0.47 23.2-22.9 

0.86x0.47 23.7-23.2 

0.96x0.47 24.2-23.7 

0.77x0.60 24.6-24.2 


Note. — ^Field of view calculated as 0(rad)=1.22/l(mm)/D(mm), where D 
is the diameter of a single antenna (12 m) and A is the observed wavelength. 


3. RESULTS 

3.1. Spatial distribution of Si-bearing molecules 

In Figs.[T]and|2l we show maps of the emission of the lines 
SiCz Jk,,k=^\6,6 - 106,5, ^®SiO 7=6 - 5 and ^^SiS 7=15 - 14 
in their ground vibrational state at different offset velocities 
with respect to the systemic velocity of the source. 

SiC 2 (Fig.[T]) displays a central component elongated in the 
NE-SW direction with a size of ~4-5" along the nebular axis 
and ~3-4" in the perpendicular direction. The elongation is 
also observed in the ^ ^SiO emission (see below) and in the 
SiO and SiC 2 maps bv iFonfrfa et al.l (120141) where the authors 
invoke a possible bipolar outflow to explain it. At the sys¬ 
temic velocity of the source and at -b6 km s“^ offset, a ring¬ 
like, clumpy and weak component is seen at ~ 10-11" from 
the central star. The angular distance between the position 
of the star and the ring stru cture, considering a distan ce of 
~ 130 pc to the star from us (iGroenewegen et all 120121) . cor¬ 
responds to a linear distance of ~2xl0^^cm. This ringlike 
component is consistent with the peak abunda nce of SiC 2 
in the outer envelope of IRC-hi 021 6 reported inlLucas et al.l 
(Il995h and the chemical model of ICernicharo et al.] (120101) . 
Although, this ringlike structure is probably Altered in our 
data given the shortest baselines used. Between the central 
and the ri nglike structure, em ission of SiC 2 is either very low 
or absent (iLucas et al.lll995h . Finally, the redshifted emission 


(-hl2km s“^) displays a quasi-circular distribution with a di¬ 
ameter of ~5-6". These brightness distributions could be in¬ 
terpreted as the SiC 2 is formed in regions close to the star, then 
it condenses onto the dust grains, and eventually it reappears 
at the outer shells of the CSE, perhaps as a hollow shell, as 
the consequence of the inter action between the UV Galactic 
radiat ion field and the CSE (iLucas et al.lll995l: lEonfrfa et~^ 

IMTl) . 

Eor ^^SiO (Eigs.[T] and [S]), the bulk of the emission 
arises from a compact central component with a size of 
2". This line also displays an extended and clumpy dis¬ 
tribution, elongated in the NE-SW and with a size of ~6- 
7". At velo cities close to the term inal expansion velocity, 
-14.5 km s~^ (ICernicharo et aO200Ql) . the brightness distribu¬ 
tion is elongated in the NE-SW direction with a size of ~3- 
4". The blueshifted emission at -13kms“^ displays a de¬ 
crease just in front of the star, which can be interpreted as self¬ 
absorption and probably absorption of the continuum emis¬ 
sion mostly coming from the star. There is no conclusive 
explanation for the elongation; nevertheless, some authors 
pointed out that it c ould evidence the p resence of a bipolar 
outflow in the CSE (lEonfria et al .1120141) . The possible pres¬ 
ence of a binary compani on to CWLeo could also play a deci¬ 
sive r ole in this scenario (iGuelin et al.ll 19931: ICernicharo et al.l 

l2ml . 

The observed brightness distributions of the vibrationally 
excited SiS lines are expected to be compact and centered on 
the star since the involved levels are excited at the high tem¬ 
peratures prevailing close to the star. The SiS /=15 - 14 v=0 
line, which displays maser emission (lEonfria Exposito et al.l 
120061) . shows a circular brightness distribution with a diame¬ 
ter of ~2" at the systemic velocity of the source. Eor the SiS 
/=15 - 14 v>l lines, the observed distributions are not spa¬ 
tially resolved. Eig.[T] shows the emission of the ^^SiS /=15 
- 14 v=0 line, which displays a brightness distribution of a 
compact source surrounding the central star with a diameter 
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of ~2". A large-scale artificial modulation can be seen below 
the 25(T level. Since it is related to the visibilities at short 
baselines, we do not expect it to modify the shape and flux of 
the compact central component of the brightness distribution. 
Also, the quality of the several SiS and ^^SiS maps is affected 
owing to the low uv coverage and non-neglible contribution 
of the sidelobes for the setups 3-5 (see Table[T]). The lines that 
lie in the range covered by the setup 6, are those of high vibra- 
tionally excited states (e.g. SiS /=15 - 14 v>10, ^^SiS /=15 

- 14 v>6) which are tentative and spatially unresolved. 

3.2. Vibrationally excited SiS 

The Cycle 0 observations with ALMA allowed us to detect 
several /=15 - 14 lines of high vibrationally excited states of 
SiS isotopologues, in particular, up to y=7 for the main iso- 
topologue (see Fig.[3]). The SiS /=15 - 14 v=8 and v=10 lines 
are probably blended with unidentified lines, so we consider 
them to be tentative. Additionally, the SiS /=15 - 14 \=9 
line is considered tentatively detected because even though its 
FWHM measured seems to follow the trend shown in Fig.[3l 
its measured integrated intensity is underestimated consider¬ 
ing the population diagram of Fig.[3l All these lines display a 
compact unresolved emission peaking at the central star. The 
FWHM of the lines measured from the spectra at the stel¬ 
lar position decreases with increasing upper level energy (see 
Fig.[3]). We verified this behavior for SiS, ^^SiS, ^^SiS, Si^^S 
and Si^^S. In the dust formation region, the gas displays a ve¬ 
locity gradient as a function of the radial distance to the star, 
i.e., the closer to the star the lower the expansion velocity 
(lAgundez et aklllM^ and references therein). Hence, those 
lines involving higher vibrational states, which are excited in 
inner and warmer regions, are narrower. Thermal broaden¬ 
ing for the SiS lines excited in the dust formation region is 
~1 km s“\ and thus this mechanism could only account par¬ 
tially for the FWHM variation of the lines. 

We analyzed the excitati on conditions of SiS with t he ro¬ 
tational diagram technique (iGoldsmith & Lang^l 19991) using 
the spectra at the stellar position (see Fig. [3]). We considered 
two different linear trends for the observational data: one for 
the transitions with Eup<2500K (i.e. v=0, 1 and 2), and a 
different trend for the rest (i.e. v=3-7). The v>8 lines are ex¬ 
cluded from the fit. The values derived from the v=0 to v=2 
fit are uncertain owing to the maser nature of the SiS /=15 

- 14 v=0 line and also due to possible optically thick emis¬ 
sion. Both data series show a linear behavior consistent with 
a single vibrational (excitation) temperature for each of the 
fits. The vibrational temperature derived from the fit involv¬ 
ing the transitions in high-energy vibrational states is higher 
than the temperature derived for the low vibrational levels. 
Therefore, the emission produced by SiS transitions in high- 
energy vibrational states can only arise from regions close to 
the photosphere. This r esult is similar to the one obtained by 
ICernicharo et all (1201 Ih for HCN. 

We used a larg e velocity gradien t (LVG) code to model 
the SiS emission (ICernicharol l2Q 1 2b . Further details about 
t he spe ctroscopic data used in the calculations are given in 
§13.2.11 The SiS collisional data were taken from iTobola et al.l 
(l2008l) and extrapolated to high rovibrational levels. We 
adopted a distance to the star of -130 pc, an effective tem¬ 
perature of -2330 K and a stellar radius of -4x10^^ cm as 
input for our model (ICernicharo et al.l |200Q|: iMonnier et al.l 
120001: iGroenewegen et al.ll2012b. We used t wo different mod¬ 
els: (0 the model of lAgundez et al.l (120121) . which was used 
to reproduce the molecular abundances in the inner layers 


of IRC-rl0216, that we call the “2012 model”, and, (//) this 
work, which is a modification of the model (/). We modified 
the H 2 density, decreas ing it by a factor -2, as described in 
ICernicharo et al.l (120131) (this was used to reproduce the dust 
nucleation zone, 1-lOR*, of IRC-i-10216). Additionally, the 
SiS abundance in the dust nucleation zone needs to be bal¬ 
anced with a similar increase (factor of -2) to avoid an un¬ 
derestimation of the SiS emission. From these models (see 
Fig. [3]) we obtain a good agreement with the vibrational tem¬ 
perature derived from the vibrational diagram and moderate 
discrepancies with the total column density within a factor 
of -2. These discrepancies in the column density may be 
explained by the dilution due to the size of the emitting re¬ 
gion compared to the half power beam width of the synthetic 
beam, which would increase the optical depth of the lines. 
The size of the emitting region should decrease with the vi¬ 
brational state owing to the energies needed to excite those 
lines. With our models, we found moderate to high optical 
depths t(v=1)-10 to t(v=4)-0.8 for the abundance profile (/) 
and t(v=1)-47 to t(v=8)-1.0 for the abundance profile (ii). 

3.2.1. SiS potential energy function and vibrational dipole moment 

The potential energy function used to describe the in- 
ternuclear motions of the SiS isotopologues is a Born- 
Oppenheimer (BO) potential properly extended to accom¬ 
modate Born-Oppenheimer breakdo wn (BOB) corrections . 
The effective po t ential has the forni (ICampbell et al. I [19931: 

Ram et al~1ll997l:lDulick et aOl998l:ICoxon & Hajigeorgiou I 

2nnnli 
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where Ma and Mb are the silicon and sulfur atomic masses 
and p is the reduced mass. The BO potential is given by 


where 


and 


y®°(r) = De 


1 - 1 ^ 
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( 2 ) 

(3) 

(4) 

(5) 


and the BOB potential and centrifugal correction terms are 
represented by the power expansions 


VA(r) = Y,4z‘ 

(6) 

Vsir) = ^ wf z' 

i=l 

(7) 


(8) 


i=l 


i=l 


The potential parameters were obtained by nonlinear least 
squares fitting to the observed infrared and microwave line 
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Fig. 3. — (Top) Central pixel spectra of the SiS 7=15 - 14 v=l, 3, 5 and 7 lines. Unidentified lines are marked with an U. Baselines are indicated in each box 
with a dashed line. (Middle) FWHM of the v>l detected lines of SiS isotopologues in the central pixel versus the upper energy of the transition. The tentative 
lines are marked with a (T). We adopted a 25% value of the total measure for the uncertainties, which takes into account uncertainties of the calibration, statistical 
errors, and the contribution of the sidelobes of the dirty beam. (Bottom) Vibrational diagram of the SiS 7=15-14 lines. The vibrational quantum number of each 
line is shown. The errorbars represent a 20% of the integrated intensity of the line and a 50% for the v=8, 9 and 10 lines. The values derived for the vibrational 
temperature and column density are shown with their formal errors for two fits: v=0-2 and v=3-7. The LVG models (see § 13.2> are plotted with a rhombus and a 
cross. 


positions of the SiS isotopologues up to v=12. The final 
data set included a total of 2863 lines, 414 rotational tran¬ 
sition s (iTieman et al. II1972I: ISanz et al. II2003I: iMuller et al. I 
2007h and 2449 r ovibrational transitions (iBirk et al. 1119721 
Frum et anil99Ql) . Mass-inde pendent Dunh a m coe fficients, 
Uij, have been derived by iMuller et al. I (l2007l) . The 
fit to the potential energy function w as performed us¬ 
ing the Levenber g-Marquardt algorithm (iLevenberg 1 11944 
iMarquardFI 1 1 963h to minimize the;^^ function, with the line 
positions weighted by the square of the experimental un¬ 
certainties. The rovibrational energy levels of the isotopo¬ 
logues needed to calculate the line positions were computed 
by solving the ra dial Schrodinger eq uation using the varia¬ 
tional method of iHarris et alTI (119651) along with harmonic- 
type basis functions. The potential parameters obtained in 
the fit are given in the comments of the Tabled The final 


value was 1.767. The unweighted standard deviations for 
the rotational and, vibrational-rotational line positions were 
0.0229 MHz and 0.000657 cm“\ respectively. 

The dipole mom ent function used for SiS was determined 
semiempirically bv iPineiro et al. I (Il987h . and the dipole mo¬ 
ment matrix elements were computed up to v=4. In Table[2] 
we provide them up to v=12. The agreemen t between our 
calculations and those of iPineiro et al. I (Il987h for v<4 is ex¬ 
cellent. 

3.3. Other Si-bearing species 

SiC was detected in IRC-f 10216 bv ICernicharo et al.l (119891) 
with line profiles indicating that the molecule was produced 
in an external shell, probably as a product of the photodis¬ 
sociation of SiC 2 . Our observations did not cover any fre¬ 
quency range where SiC lines could arise; however, lines of 
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^^SiC and SiC v=l lie in this range, but they were not de¬ 
tected. For a temperature in the photosphere of 2300 K a sig¬ 
nificant number of SiC molecules could be in the v=l state 
(Eup~1400K) and higher vibrational levels. We derived an 
upper limit to the SiC column density of 4.4x10^^ cm“^ from 
the ^^SiC upper limit, where we used an isotopic ^^Si/^^Si ra¬ 
tio of 20 (ICernicharo et al.|[T98^ 1199 ll) . Hence, SiC 2 is the 
main carrier of SiC bonds in the gas phase in the dust forma¬ 
tion zone of IRC-i-10216. 

4. CONCLUSION 

ALMA has proved to be an outstanding tool to study the 
molecular emission from CSEs of evolved stars, even at the 
early stages of its development. In particular, ALMA allowed 
us to detect SiS rotational lines in high-energy vibrational 
states that have been analyzed to constrain the physical con¬ 
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ditions of the innermost shells of IRC-f 10216. We found that 
these lines should be excited in regions close to the photo¬ 
sphere of IRC-f 10216. It also has served to unveil the dif¬ 
ferent brightness distributions of Si-bearing molecules. We 
expect that future ALMA science, with its full suite of capa¬ 
bilities ready for the next observation cycle, would give us the 
chance to map the brightness distributions of these Si-bearing 
molecules in greater detail, allowing us to understand their 
formation mechanisms. 
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TABLE 2 

Calculated vibrational dipole moment matrix elements (D) for 

v'/v” 0 _1_2_3_4_5_6_7_8_9_H3_ 11 12 

0 1.7420 

1 -0.1348 1.7658 

2 -0.6333x10-2 -0.1907 1.7896 

3 0.4068x10-3 0.1101x10-' 0.2336 1.8135 

4 -0.3379x10-^ -0.8184x10-3 -0.1564x10-' 0.2697 1.8375 

5 -0.3359x10-3 -0.7613x10-^ -0.1302x10-2 0.2027x10-' -0.3016 1.8615 

6 0.3804x10-® -0.8306x10-3 0.1329x10-3 -0.1852x10-2 -0.2493x10-' 0.3304 1.8857 

7 0.4752x10-2 0.1018x10-3 0.1569x10-'' -0.2046x10-3 0.2465x10-2 o.2961xlO-' -0.3568 1.9100 

8 -0.6394x10-" -0.1363x10-® -0.2061x10-3 0.2587x10-'* -0.2916x10-3 -0.3137x10-2 0.3433x10-' 0.3814 1.9343 

9 -0.9087x10-'* -0.1952x10-2 -0.2934x10-® 0.3612x10-® -0.3918x10-'* -0.3943x10-3 0.3866x10-2 0.3909x10-* -0.4045 1.9587 

10 0.1339x10-** 0.2936x10-" 0.4441x10-2 -0.5435x10-® 0.5779x10-® 0.5594x10-'* -0.5130x10-3 .0.4649x10-2 0.4388x10-' 0.4262 1.9833 

11 0.2005x10-'® 0.4562x10-** 0.7033x10-" -0.8654x10-2 0.9143x10-® 0.8674x10-® -0.7649x10-'* -0.6482x10-" 0.5485x10-2 0.4871x10-* -0.4469 2.0079 

12 -0.2954x10-" -0.7182x10-'®-0.1147x10-" 0.1437x10-2 -0.1525x10-®-0.1437x10-® 0.1241x10-'* O.lOllxlQ-" -0.8003x10-" -0.6373x10-2 0.5358x10-' 0.4667 2.0326 

Note. — Intemuclear potential energy parameters for the SiS molecule: D^=51993.7274cm-*, re=1.929260274 A, ;So=5.9671641, /?i =5.8793706, /?2=9.3200373, /?3=18.402479, /?4=18.451287, 
749.1309cm-* amu, u|‘=6844.274cm-* amu, Uj=-1003.321 cm-* amu, u| =7128.960cm-* amu. 



